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Abstract. We study the response functions (chromo-electric susceptibilities) for an interacting quark-gluon 
plasma. The interaction effects have been encoded in the effective fugacities for quasi-partons which are 
extracted self-consistently from the two equations of state for hot QCD. The first one is the fully pertur- 
bative 0{g^) EOS and, the second one which is 0{g^ ln{l/g)), incorporates some non-perturbative effects. 
We find that response function shows large deviations from the ideal behavior. We further determine the 
temperature dependence of the Debye mass by fixing the effective coupling constant which appears in 
the transport equation. We show that our formalism naturally yields the leading order HTL expression 
for the Debye mass if we employ the ideal EOS. Employing the Debye mass, we estimate the dissociation 
temperatures for various charmonium and bottomonium bound states. These results are consistent with 
the current theoretical studies. 
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1 Introduction 

It is expected that at high temperatures (T ~ 150 — 
200MeV) and high densities (p ~ lOGeV/ fm^) nuclear 
matter undergoes a deconfinemcnt transition to the quark- 
gluonic phase. This phase is under intense investigation in 
heavy ion collisions, and already, interesting results have 
been reported by Relativistic Heavy Ion Collider(RHIC) 
experiments [T]. As an important development, flow measure: 
suggest that close to the transition temperature Tc, the 
quark-gluon plasma (QGP) phase is strongly interacting 
— showing an almost perfect liquid behavior, with very 
low viscosity to entropy ratio — rather than showing a 
behavior close to that of an ideal gas. See ref. [3] for a 
comprehensive review of experimental observations from 
RHIC, and ref. [Tl[4l[5l|6l[7j for other recent experimen- 
tal results. On the other hand, lattice computations [H 
[S] also suggest that QGP is strongly interacting even at 
T = 2Tc- This finding has been reproduced by a number of 
other theoretical studies — by employing AdS/CFT corre- 
spondence in the strongly interacting regime of QCD|10]. 
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by molecular dynamical simulations for classical strongly 
coupled systems [TT], and by model calculations with Au- 
Au data from RHIC [T2lfT3] . In the backdrop of the above 
developments, a number of standard diagnostics, such as 
J / tjj suppression and strangeness enhancement, which have 
been proposed to probe QGP also need to be re-examined. 
It is also of importance to address other transport proper- 
ties, production and equilibration dynamics, and the phys- 
ical manifestations of pre-equilibrium evolution. 

If this is the case, as it indeed appears to be, then the 
plasma interactions would be largely in the non-perturbative 
regime; in this regime, few analytic techniques are avail- 
able for a robust theoretical analysis. Effective interac- 
tion approaches are needed. In this direction, considerable 
as already been done and we refer the reader to ref. 
[T4Ul5yl6[fT7l[T8l[T9l[20ll2T] for some of the theoretical re- 
sults. 

The effective approaches emphasize the collective ori- 
gin of the plasma properties which can be best under- 
stood within a semi-classical framework. Indeed, in a re- 
cent work [22], the successes of hydrodynamics in inter- 
preting and understanding the experimental observations 
from RHIC has been reviewed. Since more exciting and 
discerning data are expected from LHC experiments soon, 
and given the above context, it is worthwhile exploring 
semi-classical techniques to understand the properties of 
QGP in heavy ion collisions. In this context, it is known by 
now that a classical behavior emerges natu- 

rally when one considers hard thermal loop(HTL) contri- 
butions. A local formulation of HTL effective action has 
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been obtained by Blaizot and lancu who have succeeded in 
rewriting the HTL cfFcctivc theory as a kinetic theory with 
a Vlasov term [27.28, 29.30] . A significant development in 
this direction is the realization that the HTL effects are, 
in fact, essentially classical and that they are much eas- 
ier to handle within the frame work of classical transport 
equations |3 H I32 | l33]. Thus, the semi-classical techniques 
appear to hold the promise of providing tools to under- 
stand the bulk properties of QGP. 

The present paper continues the theme, and its central 
aim is to combine the kinetic equation approach which 
yields the transport properties, with the hot QCD equa- 
tions of state to make predictions which can perhaps be 
tested in heavy ion collisions. Recently, Ranjan and Rav- 
ishankar have developed a systematic approach to deter- 
mine fully the response functions of QGP, with a special 
emphasis on the color charge as a dynamical variable |18| . 
In parallel, Chandra, Kumar and Ravishankar have suc- 
ceeded in adapting two hot QCD EOS to make predictions 
for heavy ion collisions [341135] . They have shown that the 
interaction effects which modify the equations of state can 
be expressed by absorbing them into effective fugacities 
{zq^g) of otherwise free or weakly interacting quasi quarks 
and gluons. Since the analysis in ref. ([H]) was illustrated 
only for (the academically interesting) case of ideal quarks 
and gluons, it is but natural to bring the two studies to- 
gether and explore what the hot QCD EOS have to predict 
for heavy ion collisions. We take up this program in this 
paper. 

We determine the Dcbye mass as a function of tem- 
perature first by combining the quasi-particle description 
of improved pQCD EOS with the semi-classical transport 
theory formulation of response functions for QGP. Em- 
ploying the Debye mass so determined, we have estimated 
the dissociation temperatures for various quarkonia states. 
For our purpose, we consider two specific hot QCD equa- 
tions of state: The first, which we call EOSl is pertur- 
bative, with contributions up to Ofo^) [551157] . The sec- 
ond EOS has a free parameter (5, and is evaluated upto 
0[g^\og{l/g)]^. We denote it by E0S(5. 5 may be fine 
tuned to get a reasonably good agreement [35] with the 
lattice results [39], which we exploit here. Both the EOS 
are expected to be valid for T > 2Tc [38] . 

The paper is organized as follows: In section 2, we in- 
troduce the two hot QCD equations of state and outline 
the recently developed method [34] to adapt them for mak- 
ing definite predictions for QGP at RHIC and the forth- 
coming experiments at LHC. In section 3, we obtain the 
expressions for the response functions of interacting QGP 
and in section 4, we study their temperature dependence 
in detail. Moreover, we study the Debye screening and the 
dissociation phenomenon of heavy quarkonia states in hot 
QCD medium. In doing so we also relate the phenomcno- 
logical charge that occurs in the transport equation to lat- 
tice and experimental observables. We conclude the paper 
in section 5. 



2 Hot QCD equations of state and their 
quasi-particle description 

There are various equations of state proposed for QGP at 
RHIC. These include non-perturbative lattice EOS [39], 
hard thermal loop(HTL) resumed E0S[1D] and perturba- 
tive hot QCD equations of state [37ll36ll38] . In the present 
paper, we seek to determine the chromo-electric response 
functions for QGP by employing two EOS: (i) the fully 
perturbative 0{g^) hot QCD EOS proposed by Arnold 
and Zhai[36j and Zhai and Kastening [37], and (ii) The 
EOS of 0[g^{\n{l/g) + 6)] determined by Kajantie et al 
[38j , by incorporating contributions from non-perturbative 
scales, gT and g^T. We employ the method recently for- 
mulated by Ranjan and Ravishankar [18] to extract the 
chromo-electric permittivities of the medium. The EOS 
which we label E0Si5 is given by 
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EOSl is obtained from this equation by dropping the last 
term which has contributions of 0{g^{ln{l/g) +5)). The 
phenomenological parameter d is introduced in [38j to in- 
corporate the undetermined contributions oiO{g^). It also 
acts as a fitting parameter to get the best agreement with 
the lattice results. . 

EOSl and E0S(5 have several ambiguities, associated 
with the renormalization scale (/^ms), the scale parameter 
yl-r/^Ms which occurs in the expression for the running 
coupling constant a^^ and the value of the phenomeno- 
logical parameter S. The arbitrariness in fixing (nffg) has 
been discussed well in literature and a popular way out is 
the BLM criterion due to Brodsky, Lepage and Macken- 
zie [H]. In this criterion, the value of {hms) allowed to 
vary between ttT and AttT [42]. Here, we choose Hfjs ~ 
2.157rT w 6.752T[ig close to the central value 2ttT. With 
this particular choice, all the contributions due to the 
logarithms containing /ij^/g are very small. For the scale 
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Fig. 1. (Color online) Relative equation of state wrt ideal EOS 
for pure gauge theory plasma as a function of T/Tc for various 
values of 5. 




T 



Fig. 2. (Color online) Relative equation of state wrt ideal EOS 
for full QCD plasma with Nf = 2, 3 as a function of T /Tc for 
various values of 5. 

parameter At, we follow Huang and Lissia [H] and set 
At/Aj^s = exp(7ij + l/22)/47r k. 0.148. They found that 
this choice makes the coupling g^{T) is optimal for lattice 
perturbative calculations. The same value has also been 
employed by others, see e.g. [381143]. 



2.1 The underlying distribution functions 

The construction of the distribution functions that under- 
lie the EOS, in terms of effective quarks and gluons which 
act as quasi-excitations, has been discussed by Chandra 
et. aZ[34] in the specific context of EOSl and E0Si5. To 
review the method briefly, all the terms that represent 
interactions are collected together by recasting them as 
effective fugacities (zq^g) for the otherwise free quarks and 
gluons. Of course, the pure gauge theory case is simply 
obtained by putting the number of flavors, Nj = in the 
EOS. Thus, Zg represents the self interactions of the glu- 
ons, while Zq encapsulates the quark-quark and the quark 
gluon interaction terms. The quantities Zg/q can be de- 
termined from the two EOS self-consistently. In this pro- 
cedure, all the temperature effects are contained in the 
effective fugacities z = z{as(T/Tc)), where we display the 
dependence on the temperature and coupling constant ex- 
plicitly. It has been shown in ref. [34] (where the details 



can be found) that one can trade off the dependence of the 
effective fugacities on the renormalization scale ( ^J■Ms) ^y 
their dependence on the critical temperature Tc- For that 
purpose, one utilizes the one loop expression of as(T) at fi- 
nite temperature given by [44] . It should be borne in mind 
that the effective fugacity Zg/q introduced here has noth- 
ing to do with the usual fugacity which corresponds to the 
conservation of particle number. It is, therefore, unrelated 
to the baryon chemical potential in the case of quarks. 
Zg/q merely encodes the interaction effects present in hot 
QCD EOS. Its deviation from unity signifies the presence 
of non-ideal terms in the EOS. 

In a recent work Chandra and Ravishankar [35] showed 
that the effective fugacities Zg/g modifies the dispersion 
relations. The modification captures the non-zero trace 
anomaly effects in hot QCD. They further employed this 
model to determine the shear viscosity and shear viscosity 
to entropy ratio as a function of temperature. This model 
has further been generalized to the pure lattice gauge the- 
ory EOS in [45] and the authors found that the model 
works remarkably well for lattice QCD. This quasi-particle 
description is analogous to Landau's theory of Fermi liq- 
uids. The authors [15] have further determined the tem- 
perature dependence of gluon quenching parameter and 
studied the shear viscosity (77), and its ratio with entropy 
(77/5) employing the recent work by Asakawa, Miiller and 
Bass [37] and Majumder, Miiller and Wang [35]. 

In figs. 1 and 2, we display the behavior of E0S(5 for 
various values of the parameter 5. The figures show the 
pure gauge theory contributions to the EOS and full QCD 
separately. We remark parenthetically that the studies in 
the earlier work [M] were confined to EOSl and the spe- 
cial case (5 = in E0S(5. For the details on EOSl and 
E0S(5 for (5 = 0, we refer the reader to ref. [33] (see fig.l- 
7 of ref. [33]). First of all, we see that as S increases in 
magnitude, the EOS, for both pure gauge theory and full 
QCD, become softer, with P/Pi taking smaller values, we 
denote the the ratio P/Pj by Ri. Kajantie[35| obtains 
the best fit with the lattice results of Boyd et. al.[32| by 
choosing a value S = 0.7. We find that to get agreement 
with the more recent results of Karsch [39], 5 « 1.0 is 
preferred, when we consider T > 2Tc. Over all, we find 
that the range of values 0.8 < 5 < 1.2 gives a reasonably 
good qualitative agreement with the lattice results for the 
screening lengths. Finally, we set Aj^jg = Tc, which is close 
to the value 0.87rc found by Gupta [5D]. Here, we wish to 
mention that there is an uncertainty in fixing the free pa- 
rameter S. This follows from the freedom in choosing the 
QCD renormalization scale at high temperature. This has 
been investigated in detail by Blaizot, lancu and Rcbhan 
[51] . The value of i5 in the present paper has been obtained 
by employing the one loop expression for the running cou- 
pling constant and the QCD renormalization scale deter- 
mined in ref. [44] . 

The behavior of the corresponding fugacities, as a func- 
tion of temperature, is shown in fig. 3. It may be seen 
that < Zg,g < 1.0 which ensures the convergence of the 
method to determine the effective fugacities from the hot 
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QCD EOS. Wc now proceed to determine the response of 
the plasma in the next section. 




2 2.5 3 3.5 4 4.5 5 



T/T, 

Fig. 3. (Color online) Effective parton fugacities {zg,q) quarks 
determined from EOSiJ as a function of temperature. Note that 
the behavior is shown for 5 = 1.0. 



3 Response functions for interacting qgp 

Recently Ranjan and Ravishankar [18| have determined 
the form of chromo-electric response functions for collision 
less quark-gluon plasma within the framework of semi- 
classical transport theory. They have set up the transport 
equation in the extended phase space including the SU(3) 
group space corresponding to dynamical color degree of 
freedom. They have taken the distribution function in a 
coherent state basis defined over the extended single par- 
ticle phase space TZ^ (g) Cg, where Cg = Q /TL \& the phase 
space corresponding to the color degree of freedom, ob- 
tained as a coset space by factoring the group space by 
the stabilizer group Tl of any reference state in the Hilbert 
space. Having been employed to study the ideal case, the 
formalism has not been applied to examine the behavior 
of the plasma with a realistic EOS. We employ the results 
of the previous section and rectify this drawback, by in- 
corporating the interaction effects as represented by EOSl 
and E0S5. 

A brief comment on the response functions. In contrast 
to electrodynamic plasma, the chromo-electric response 
has a richer structure. Apart from the standard permittiv- 
ity which we shall call Abelian and denote by eA^ there are 
additional response functions, their number depending on 
the color carried by the partons. Thus, quarks have an ad- 
ditional response function which affects the non- Abelian 
coupling. The corresponding permittivity will be called 
non- Abelian, and denoted by e^v- The two functions ex- 
haust the response in the quark sector. The gluonic sector, 
arising from the adjoint representation of the gauge group 
admits yet another kind of response, corresponding to ten- 
sor excitations. These excitations are not allowed in the 
quark sector (which emerges from the fundamental repre- 
sentation of the gauge group). We consider each of these 



response functions for the interacting QGP. The response 
functions are obtained in the temporal gauge. 

Consider first the familiar Abelian component of the 
response tA- For an isotropic plasma(in the absence of 
chromo-magnetic fields), its expression is given by [TH] 

eA(c^,fe) = l + Q2/„(^^fc) (2) 

where = Q°^Q'^ is the color charge magnitude squared, 
and /q is determined by the equilibrium distribution func- 
tion thus: 

J OJ - ^ + IT op., 

The non-Abelian response function, which has been 
evaluated in the long wavelength limit, is given by 

+ (3) 

UJ 

where Ii is defined as 

^ ' 3 yj (cj-^) dp, ^) 

We recall that the new constitutive Yang-Mills equations, 
in the presence of the medium, are given by 

p^iu, k) + iQ'^E^iuj, k)kJo{'^, k) 

r^l £alra r 

.'±±_ / I^{J^k')\^,^^A\{,^-uo'^k-k') 
^ J 

xE^{uj',k')duj' d^k' = 0. (4) 

j;(^, k) + iQ^EIiio, k)S,, hiuj, k)\^^, 

= 0. (5) 

As pointed out in [TS], the Abelian and non-Abelian 
responses are not independent of each other. Gauge in- 
variance relates them, by virtue of which we can obtain 
both from a common generating function as follows: 

I d f k ■ p 

^° " k^d^ J ^"^^"^ —)kidpjegd^p 

= A^'H I '"^"^ ^)9pJe,d'p) ■ (6) 

We further recall that these expansions are determined 
when the system is displaced slightly from its equilibrium, 
in the coUisionless limit. 



3.1 Ideal response 

It is convenient to first write the expressions for the re- 
sponses of ideal distributions for quarks and gluons. The 
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responses due to EOSl and E0S5 get a simple modifica- 
tion over their ideal forms since we have mapped success- 
fully the interaction effects into quasi free partons with 
effective fugacities. Thus, in the ideal case we have, for 
the quarks, 



- k 



2}] 



(7) 



and the non-Abelian response function is given by 



9 



(8) 



gu{zg) has the series expansion 



Note that 52(1) — 1 gives the ideal limit. 

Similarly, the corresponding expressions in the quark 
sector are obtained as 



uj + k 



- k 



2}] 



(12) 



and the non-Abelian response for effective quarks reads: 



The imaginary part of Abelian(eyi) and non-Abelian 
component (cat) of the chromo-electric permittivity can be 
easily evaluated by the standard Landau it prescription. 
These are needed to obtain landau damping which we do 
not study here. The gluon effective charge(Qg) and the 
quark effective charges (Q^ are related to each other by 
the total number of color Nc for SU{Nc) as Qg = NcQg = 
NcQ^, where Q is the color charge of a single quark. It is 
easy to see that the gluonic susceptibility, x^''' = ^ — 1 

is identical to the quark susceptibility, x'"''^ 
they are related as, 



N 



(9) 



3.2 Interaction effects 



eN = { 



9 



UJUJ' 



7}- (13) 



The function J2{zq) = ^ f2{zq) where f2{zq) is defined 
via the integral below. 



■dx = r{v)f„{zq) 



Zq ^exp{x) + 1 

00 I 



1=1 



and /^(l) = 1. 



4 Effective charges and relative 
susceptibilities 



We now consider the modification that the above expres- 
sions undergo permittivities arising because of the new 
EOS. Recall that the corresponding equilibrium distribu- 
tion functions differ from each other only in their form for 
the chemical potentials fJ-q^g. The responses thus depend 
on the interactions implicitly through an explicit depen- 
dence on Zq^g. 

Considering the gluonic case, i. e., pure gauge theory 
first, we get the expressions for the two permittivities as 



cu + k 



k 



and the non-Abelian response function as 



, 87:^ N.Q^T^g'.jzg) 1 
= {1 n 

9 LUUJ' 



2}], (10) 



(11) 



The function g2{zg) — ■§7g2{zg) where g^izg) is defined 
via the integral below. 



lo Zg exp{x) — 1 



dx = r{iy)g^{zg). 



eq. () 1011131 admit a simple physical interpretation, when 
compared with their counterparts eq. ([7] -[5]) . Indeed, the 
sole effect of the interactions on the transport properties is 
to merely renormalize the the quark and the gluon charges 
Qg^q as shown below: 



l92{Zg); Ql 



l2 fl 



f2iZq). 



The renormalization factors g2{zg), f2{zq) further possess 
the significance of chromo-electric susceptibilities, relative 
to the ideal values. To see that, we note that the Abelian 
and the non-Abelian strengths for gluons as well as quarks 
suffer the same renormalization reflecting the underlying 
gauge invariance. Furthermore, the expressions for the rel- 
ative susceptibilities are given by. 



and 



X{z) ^ A(z)_ _ Af{z)_ 
x(i) - ^(1) ~ U{1) 



fl^izq) for quarks, 
92{zg) for gluons ^ ' 



X^'Hzq) ^ A^^Hzq) ^ U^'>){Zq) ^ Nf f^jZg 
X^oKzg) A(9){zg) M3){Zg) 2iVe^ 



■15) 



Note that the relative susceptibilities are entirely func- 
tions of the single variable T/Tc, and are independent of 
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(oj, k). The dependence of the susceptibilities on (w, k) has 
already been studied in detail in ref.[18j. We merely con- 
centrate on the temperature dependence below. 

Before we go on to discuss the susceptibilities and other 
bulk properties, we point out an essential care to be taken 
in using the above susceptibilities for determining the re- 
sponse of the plasma. For pure gauge theory, only the glu- 
onic part contributes, while for the full QCD. we have to 
necessarily take the contribution from both the quark and 
the gluonic sector. We discuss both the cases below. The 
response functions for the full QCD is obtained by sum- 
ming over the response function for quark as well as gluon 
plasma. The relative susceptibility for full QCD plasma is 
given by 



7^' 



iV//^(z,)+27Veg^(z,) 
Nf + 27Vc 
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Fig. 5. (Color online) Relative susceptibility, defined 
(|16|) . for the full QCD plasma as a function of T/Tc, for 
and EOS5 (3 = 1). We have studied the cases Nf = 2,3 



m eq. 
EOSl 



(16) 



4.1 Behavior of the susceptibilities 

We now proceed to study the behavior of the relative sus- 
ceptibilities displayed in eqs. p4|) . p5|) and p6|) as func- 
tions of temperature. As observed, relative susceptibilities 
for both quarks and gluons scale with T/Tc- We have plot- 
ted the relative susceptibilities 72,, Ti-qg and TV as functions 
of r/Tc (See figs.4-6), for both EOSl and Y.0S5. Please 
note that we have chosen 5 = 1.0 in E0Si5. 

Fig. 4 shows the relative susceptibility of a purely glu- 
onic plasma as a function of temperature for EOSl and 
E0S(5. 



of the phenomenological coupling Q^, 
value. 



relative to its ideal 



The behavior of quark gluon plasma is not qualita- 
tively different from that of a purely gluonic plasma, as 
may be seen from fig. 5. In other words, the quark contri- 
bution is of the same order as the purely gluonic contribu- 
tion. However, the relative contribution from the quarks 
and the gluons does depend on the EOS considered. In- 
deed, with EOSl (where interactions up to 0{g^) are in- 
cluded), fig. 6 shows that the gluonic contribution domi- 
nates over the quark contribution for both Nf = 2 and 
Nf = 3 and same is true for E0S(5. The dominance is 
more pronounced for EOSl as comapare to EOSJ. 

At this juncture, we clarify that by ideal we mean that 
we employ the ideal EOS . The expressions so obtained 
are the same as the one loop result, or equivalently, the 
leading order HTL result. Similarly, the phrase interacting 
QGP refers to the contribution of the non-ideal terms in 
the hot QCD EOS, which is equivalent to the inclusion of 
higher order corrections to the one loop result. 



Fig. 4. (Color online) Relative susceptibility, 32 (-^s) (see eq. 
I|14p . for pure gauge theory plasma as a function of T/Tc for 
EOSl and EOS5 ((5 = 1). 



We see From fig. 4 that the susceptibility of a purely 
gluonic plasma is weaker in the presence of interactions, 
approaching its ideal value asymptotically with increasing 
temperatures. Equivalently, there is a decrease in the value 




Fig. 6. (Color online) Ratio of the quark to gluonic contribu- 
tions to the susceptibility (see eq. (|15|l as a function of T/Tc for 
EOSl and E0S(5 for 5 = 1 
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4.2 Debye mass from the chromoelectric susceptibility 

The static or zero frequency limit of the response function 
leads to the well known relation with the Debye mass as 

follows, e(fc,r) = 1 + 

One can compute the Debye mass for full QCD by 
utilizing the additive nature of Abelian response functions 
of partons in QCD. The full permittivity for hot QCD 
plasma will be obtained, 



(17) 



which naturally incorporates the additivity of Debye mass 
square in thermal QCD. 

For an ideal QGP the Debye mass comes out to be. o.25 



and fig. 8 respectively. To compare these results with lat- 
tice predictions of Zantow and Kazmarek[S2], we choose 
Tc = 0.27 GeV for pure gauge theory and = 0.203, 0.197 
for 2- and 3-flavor QCD respectively. The screening length 
for both EOSl and E0S2 show qualitative agreement with 
these results in [52]. As expected, the agreement is more 
pronounced for E0S5 as compare to EOSl. We shall now 
proceed to determine the dissociation temperatures for 
various quarkonia states in a hot QCD medium. 



0.35 
0.3 



(M^y = 87r^Q^T^{Nc/3 + Nf/6). 



(18) 



Interestingly, this matches with the leading order HTL 
expression for the Debye mass if one identifies the phe- 
nomenological charge Q as 



(19) 



M 



.Nr. 



N 



To see how the interactions influence the Debye screening. 



we consider the ratio, Raij^ 
given as. 



R 



2N,g',{zg) + Nff^izg 
2Nc + Nf 



0.2 



0.15 



0.1 




Fig. 7. (Color onhne) 1/Md as a function of T/Tc for EOSl 



where g{T) is the QCD running coupling constant at fi- 
nite temperature. This fixes the phenomenological charge 
Q. We employ this to study how the inclusion of interca- 
tions influence the Debye screening and the dissociation 
phenomenon of heavy quarkonia in a hot QCD medium. 

The current approach to determine the Debye mass is 
rather different from that in ref . |34j . Here, we have deter- 
mined the Debye mass from the chromo-electric permit- 
tivity recently obtained by Ranjan and Ravishankar[T5l . 
On the other hand, the method adopted in ref. [34] is based 
on the earlier works of Kelly et. a/ |311l32| . We remark that 
the two approaches are equivalent at the one loop level 
provided that the phenomenological charge is fixed ap- 
propriately. 

For the interacting QGP the Debye mass can be writ- 
ten by. 




(20) 



{MofKMjjf, which is 



(21) 



The expression of Ri\[,j is the same as that for relative 
susceptibility TV which follows from the expression for the 
Debye mass of full QCD in terms of the response function 
( egfTTj). It is clear from figs. 4 and 5 that the inclusion of 
the interactions significantly lower the Debye mass. The 
Debye mass will approach its ideal value only asymptoti- 
cally. We have shown the behavior of the screening length 
as a fimction of temperature for EOSl and E0S(5 in fig. 7 



3.5 4 

Fig. 8. (Color online) 1/Md as 'a function of T/T^ for EOS(5. 



4.3 Dissociation temperatures 

To determine the dissociation temperatures of heavy quarko- 
nia, we follow a well known criterion-whenever the rms ra- 
dius, rqq > 1/M]j for a particular heavy quarkonia bound 
state, the state will dissociate in the medium. The equality 
sign will yield the dissociation temperature. To this end, 
we shall employ the rms radius of various bound states of 
charmonium and bottomonium in ref. [53] . We have shown 
the dissociation temperatures for various quarkonia in Ta- 
ble I along with the fiavor dependence. The dissociation 
temperatures decrease with increasing number of flavors. 
The J/'Z' dissociation temperature varies between l.QTc to 
1.5Tc for E0S<S and 2.28Tc to IMT^ for EOSl. Note that 
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Table 1. The dissociation temperature(T£)) for various 
quarkonia states (in unit of Tc). Note that we have employed 
one loop expression for the running coupling constant [44] and 
Tc =0.27, 0.203 and 0.197 for Nf = 0, Nf = 2 and Nf ^ 3 
flavor QCD [52 respectively. 



Hot EOS 


Quarkonium 


Nf=0 


Nf^2 


Nf=3 


EOS(5 




1.92 


1.60 


1.50 






< 1.50 


< 1.50 


< 1.50 




T 


3.88 


3.30 


2.77 






1.70 


1.50 


< 1.50 




Xc 


< 1.50 


< 1.50 


< 1.50 




Xb 


2.21 


1.87 


1.60 


EOSl 


J/w 

^' 

r 

Xb 


2.28 

<1.50 

4.35 

2.05 

< 1.50 

2.60 


1.85 

< 1.50 
3.74 
1.64 

< 1.50 
2.15 


1.58 

< 1.50 
3.22 
1.50 

< 1.50 
1.84 



the leading order HTL result for an ideal QGP. The viabil- 
ity of the two EOS, especially E0S(5 is thus phenomeno- 
logically well supported. It would also be of interest to 
extend the analysis to other signatures like strangeness 
enhancement, and also for QGP with a finite baryonic 
chemical potential [58l[59] and the HTL and HDL equa- 
tions of state|5H[60]. 
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